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Abstract 


This paper is in continuation of the authors’ recently published paper (Journal of Mathemat¬ 
ical Physics 55(2014)083519) in which computational solutions of an unified reaction-diffusion 
equation of distributed order associated with Caputo derivatives as the time-derivative and 
Riesz-Feller derivative as space derivative is derived. In the present paper, computable solutions 
of distributed order fractional reaction-diffusion equations associated with generalized Riemann- 
Liouville derivatives of fractional orders as the time-derivative and Riesz-Feller fractional deriva¬ 
tive as the space derivative are investigated. The solutions of the fractional reaction-diffusion 
equations of fractional orders are obtained in this paper. The method followed in deriving the 
solutions is that of joint Laplace and Fourier transforms. The solutions obtained are in a closed 
and computable form in terms of the familiar generalized Mittag-Leffier functions. They provide 
elegant extensions of the results given in the literature. 
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1. Introduction 


Distributed order fractional reaction-diffusion systems are studied, among others, by Saxton 
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[54,55,56], Haken [16], Langlands [28], Sokolov et al. [57-59], Saxena and Pagnini [51], Guo and 
Xu [15], Huang and Liu [24], and recent monographs on the subject [4,22,26,27,31,38]. General 
models for reaction-diffusion systems are discussed by Wilhelmsson and Lazzaro [64], Henry and 
Wearne [19,20], Henry et al. [21], Haubold et al. [17,18], Mainardi et al. [29,30], Saxena [42], 
and Saxena et al. [43,44,45,47]. Stability in reaction-diffusion systems and nonlinear oscillation 
phenomena have been discussed by Gafiychuk et al. [12,13]. Pattern formation in reaction- 
diffusion related to physical and biological sciences can be found in the work of Murray [34], 
Gross and Hohenberg [3], and Nicolis and Prigogine [36]. Recently, Engler [7] discussed the speed 
of spread for fractional reaction-diffusion. Distributed order sub-diffusion is discussed by Naber 
[35]. General models for reaction-diffusion systems are discussed by Hilfer [22], Wilhelmsson 
and Lazzaro [64], Henry and Wearne [19,20], Henry et al. [21], Mainardi et al. [29,30], Haubold 
et al. [17,18], Diethelm [4], Kuramoto [27], Hundsdorfer and Verwer [25], and Saxena et al. [42- 
48]. The fundamental and numerical solution of a reaction-diffusion equation associated with 
the Riesz fractional derivative as the space derivative is derived by Chen et al. [2]. Reaction- 
diffusion models associated with Riemann-Liouville fractional derivative as the time derivative 
and Riesz-Feller derivative as the space derivative are recently discussed by Haubold et al. [18]. 
Such equations in case of Caputo fractional derivative are recently solved by Saxena et al. [50]. 

In connection with the evolution equation for the probabilistic generalization of Voigt profile 
function, it is shown by Pagnini and Mainardi [37] that the solution of the following integro- 
differential equation 


— = oDo^^^Nix, t) + t),N{x, 0) = (1.1) 

is obtained in terms of the Fourier transform, where oDo°‘^ and are the Riesz fractional 

derivatives of orders ai and a 2 respectively, and S(x) is the Dirac-delta function, which is given 
in [37, p.l593]. Consider the Fourier transform with parameter k: 

N*{k,T) = exp[—r(|fc|“^ -I- |fc]“^)],r > 0. (1.2) 

This has motivated the authors to investigate the solutions of partial differential equations (2.1), 
(4.1) and (5.1), listed later on. The solutions are obtained in a closed and computable form in 
terms of the familiar Mittag-Lefller functions. Some known and unknown results associated with 
fractional reaction-diffusion equations and fractional reaction-diffusion of fractional orders can 
also be derived, as special cases of our findings. It may be observed that in case of distributed 
order fractional reaction-diffusion equation, the solution can be written in a compact and closed 
form in terms of a generalization of Kampe de Feriet hypergeometric series in two variables. Due 
to general character of the derived results, the results given earlier by Chen et al. [2], Haubold 
et al. [17] and Pagnini and Mainardi [37], Saxena et al. [48, 50], and others, readily follow as 
special cases of our investigations. The solutions are obtained in forms suitable for numerical 
computations. 

2. Solution of the Fractional Reaction-Diffusion Equation 

In this section, we will derive a computable solution of the one-dimensional fractional 
reaction-diffusion equation, given below in (2.1), containing generalized Riemann-Liouville frac¬ 
tional derivatives as the time-derivatives and a Riesz-Feller fractional derivative as the space 
derivative. The results obtained are in a compact and computable form in terms of the gener¬ 
alized Mittag-Lefller function, defined by (A3) in the form of the following theorem. 
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Theorem 1. Consider the one-dimensional fractional reaction-diffusion equation of fractional 
order 


t) + t) = r]^De°‘N{x, t) - ujN{x, t) + U{x, t), (2.1) 

Here uj,r],t > 0,x € R; a, 0, 71,72, ^1, <^2 are real parameters with the constraints 

1 < 71 < 2,0 < ( 5 i < 1; 1 < 72 < 2,0 < ^2 < 1, 0 < a < 2; (2.2) 

are the generalized Riemann-Liouville fractional derivative operators defined 
by (A9) with the conditions 

= fi{x); = 51 (x) 

= /2(x); ^ 

lim N{x,t) = 0. (2-3) 

|ai|—^00 

Further, uj is a constant with reaction term, xDe°‘ is the Riesz-Feller fractional derivative of 
order a and symmetry 9 defined by (All) with \6\ < min(a, 2 — a), 77 is a diffusion constant and 
U[x, t) is a nonlinear function belonging to the area of reaction-diffusion. Then the solution of 
(2.1), under the above conditions, is given by 


°° (—nY 

N{x,t) = t 7 i+- 5 i( 2 - 7 i)- 2 y^ 


_|_ ^7 i +<^ i (2-7 i )-1 7(71-72)1’ 


r—0 


27T 


X g( (k) )dfc 

°° (—nY 1'°° 

r=0 

^ ’^7i.7i + (7i-72)i'+<52(2 - 72)-1^ 

_^q77i+-52(2 - 72)-1 ^ (Z^ f g*(J,-^ 

°° (—nY F 

/ ^7i+(')'i-72)r-i / u*{k,t-OeM-ikx) 
ZTL 27r Jo J -00 


r =0 

^r +1 


27r ./o 

where )R( 7 i) > 0 , 5 R( 72 ) > 0 and )R( 7 i — 72 ) > 0,6 = w + pipai^)- 


(2.4) 


Proof: If we apply the Laplace transform with respect to the time variable and the Fourier 
transform with respect to the space variable x, use the initial conditions and the formula (All) 
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and (A12), then the given equation transforms into the form 

+ as'>'=iV*(fc,s) 

= T]'iPi{k)N*{k, s) - ojN*ik, s) + U*{k, s) 


where according to the convention followed, the symbol ^ will stand for the Laplace transform 
with respect to the time variable t and * represents the Fourier transform with respect to the 
space variable x. Solving for N*{k, s) it gives 


+ fUk) + + U*is) 


(2.5) 


where b = uj + r]tp^{k). Inverting the Laplace transform by using the result (Appendix A(19)) 
we obtain 


N*{k,t) = /*(fc)i7i+5i(2-7i)-2^('_^)r^(7l-72)r 

r—0 

^ ‘^7i>7i + (7i-72)»’+'5i(2-7i)-1^ ^ 

oo 

+ gj((fc)tT'i+‘^i(2-7l)-l ^(_o)^i(7l-72)r 

r—0 

^ '^7i>7i + (7i-72)»’+(5i(2-7i)( ) 

OO 

r—0 

^ "^7i>7i + (7i-72)»’+(52(2 - 72)-i( ) 

OO 

+ ag2(fc)ri+‘5=^(2-72)-l ^(_Q)^i(7l-72)r 

r—0 

^ '^7i>7i + (7i-72)»’+(52(2 - 72)( ^ 

pt oo 

•^0 r=0 

X <;. + (7.-72).(-^^"^)- 


l- 72 )»’-l 


( 2 . 6 ) 


The required solution to (2.1) is now obtained by taking the inverse Fourier transform of (2.6). 


3. Special Cases 


When w = 0 Theorem 1 reduces to the following: Corollary 1.1. Consider the one¬ 
dimensional fractional reaction-diffusion equation of fractional order 

Df^’^^N{x, t) + aDf'^^N{x, t) = r]o,D 0 ‘^N{x, t) + U{x, t), (3.1) 

for r],t > 0,x € R;a, 9, 71 , 72 , i^i, <^2 are real parameters with the constraints 

1 < 71 < 2,0 < (5i < 1; 1 < 72 < 2, 0 < (52 < 1,0 < a < 2; (3.2) 
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with the conditions 


^ ^^(1-5O(2-7i)^(^^ 0+) = 5i(a:) 

= / 2 (x); A/a-^.)(2-7.)^(^^0) ^ 

lim N{x, t) = 0. 

|a:|—^oo 


(3.3) 


xDq^^ is the Riesz-Feller fractional derivative of order a and symmetry 9 defined by (All) with 
|0| < min(a, 2 — a), ry is a diffusion constant and U{x,f) is a nonlinear function belonging to the 
area of reaction-diffusion. Then the solution of (3.1) is given by 


N{x,t) = t 7 l+- 5 l( 2 - 7 l)- 2 ^ 


r—0 


i-aY 

27T 


fli 72)r-y*^yj,^ 


X i^;+7l +(71 -72).+5 i (2-71)-1 Yc ik))dk 

_|_ y7l+<^l(2-7l)-l ( (~Q^) f f'n-'t2)r 

-n 2,7T J — QQ 


r=0 


°° (—nY 

_^^^7i+52(2-72)-2^i^ / <(^i-'>'=)’'/2(*)exp(-ifcx) 

r=0 d-oo 

_^q^ 7 i+- 52(2 - 72)-1 ^ (z^ f g *( JY ^ 

-n —OO 


r=0 




E 

r—Q 

Tir+l 


{-ay 


27T 


^7i + (7i-72)'r-l 

f7l,/,e 


/ OO 

U*{k, t — f) exp(—ifcx) 

-OO 


X +(^,_^,)7-^ec(fc))dfcd7 

where 9 fJ( 7 i) > 0 , 5 ft( 72 ) > 0 and (R (7 — 72 ) > 0 , 6 = vYaik). 


(3.4) 


If we set 9 = 0 then by virtue of the relation (A14), the Riesz-Feller space derivative reduces 
to a Riesz derivative and Theorem 1 reduces to the following: 

Corollary 1.2. Consider the one-dimensional reaction-diffusion equation of fractional order 
DY’^^N{x, t) Y aDY’^^N{x, t) = r]xDo°N{x, t) - ujN{x, t) + U{x, t), (3.5) 

for uj,r],t > 0,x € R; a, 9, 71 , 72 , ^i, S 2 are real parameters with the constraints 


1 < 71 < 2, 0 < di < 1; 1 < 72 < 2,0 < ^2 < 1, 0 < a < 2; 


(3.6) 
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with the conditions 


^ ^^(i- 5O(2-7O^(^^0+) = 5i(a:) 

= /2(x); A/a-^.)(2-7.)^(^^0) ^ 

lim N{x,t) = Q. (3.7) 

|a:|—^oo 


Lu is a constant with the reaction term, xDo°‘ is the Riesz fractional derivative of order a defined 
by (All), rj is a diffusion constant and U{x,f) is a nonlinear function belonging to the area of 
reaction-diffusion. Then the solution of (3.5) is the following: 


°° (—nV 1'°° 

7V(a;,t) =Fi+‘5i(2-7l)-2y- I- L / 


X + I*r))dfc 

oo 

_|_ ^7l+<5l(2-7l)-l 


/ \r /*oo 


r—0 


X g({k) exp{-ikx)Ef^+^^^_^^^^^g^^.^_^^^{-r]t'^^{ix + \k\‘^))dk 
oo (_„Y 1'°° 

+ ai'^i+'5"(2-72)-2 ^ / i(7i-72)ry*(^) exp(-ifca;) 

X (a; + |fcr))dfc 

_^Qi7i+-52 ( 2 - 72)-iy^ -L / 

—' 27r / „ 

X + |fcr))dA: 

°° / Nr pt poo 

—-- / ^7i+(7i-72))'-i / U*{k,t — ^)exp{—ikx) 

_n ^TT Jg J-oo 


r =0 

pr+l 


X _,,)y-^e(^ + ifcr))dfcde, 


(3.8) 


where 9 fJ( 7 i) > 0 , 5 R( 72 ) > 0 and 5 R (7 — 72 ) > 0 . 

If we set (5i = (52 = 0 then the generalized Riemann-Liouville fractional derivatives and 

£)i 2 ,S 2 ]-g(jyQg respectively to the Riemann-Liouville fractional derivatives and 

defined by (A5) and we arrive at the following: 


Corollary 1.3. Consider the one-dimensional fractional reaction-diffusion equation of frac¬ 
tional order 


O^Dt'^^Nix, t) + a^^Dt^^N{x, t) = T]xDg^N{x, t) - u:N{x, t) + U{x, t), (3.9) 

for uj,r],t > 0,x € R; a, 9, 71 , 72 , are real parameters with the constraints 


1 < 71 < 2; 1 < 72 < 2; 


(3.10) 
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where and are the Riemann-Liouville fractional derivative operators defined by 

(A5), with the conditions 

^( 7 i- 2 )^(^^ 0 )) = ffix);D['^^~^''N{x,0+) =gi{x) 

Di-^^-^^N{x,0) = f2ix);Di^^-^^Nix,0)=gfix), 

lim N{x,t)=0. (3-11) 

|a:|—>-oo 

where ^^iV(x,0) denotes the (71 — 2)th derivative of N{x,t) evaluated at t = 0 and other 

quantities are as defined before and the conditions on the parameters, including the one on 6 
remain the same. Then the solution of (3.9), under the above conditions, is given by 


N{x,t)=t^^-^^ 


r—Q 


{-ay 

27T 


fill exp(—ifcx) 




c'-'Y. 


r =0 


27r 


^C7i-72)r 


' —00 
-ir+l 


X gl(k) 


afn-^Y 


r—0 


i-ay 

27T 


fill 




00 (_ Y r°° 
r=0 ^ 

X eyip{—ikx)E'"^^ , . (—bEfidk 

“ (—nY f* 

+ — / / U*{k,t — f)exp{—ikx) 

r=o “^0 

x<y+(,,_,2).(-^e)dMe, 


(3.12) 


where 9 fJ( 7 i) > 0 , 5 ft( 72 ) > 0 and 5 R( 7 i — 72 ) > 0 , 6 = w + 

If we set di = (52 = 1 then the generalized Riemann-Liouville fractional derivatives 
and reduced respectively to the Caputo derivatives ^and ^Dfi^ defined by (A6) and 

we arrive at the following: 


Corollary 1.4. Consider the one-dimensional reaction-diffusion equation of fractional order 
Q N{x,t) + Oq N{x,t) = gxDe°'N{x,f) — ujN{x,t) + U{x,t) (3.13) 


where uj,g,t > o, x € i?; a, 71,72 are real parameters with the constraints 1 < 71 < 2,1 < 72 < 
2, 0 < q; < 2 with 


N{x,0+) 


f{x), —N{x,0+) = g{x)) 


(3.14) 
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and uj is a constant with reaction term, is the Riesz-Feller space fractional derivative of 

order a and symmetry 0 defined by (All), with |0| < min(a, 2 — a), 77 is the diffusion constant 
and U{x,t) is a nonlinear function belonging to the area of reaction-diffusion. The solution of 
(3.13), under the above conditions, is given by the following: 

°° (—nY 

N{x,t)='S^ - / f({k) exp{—ikx) 

h'o 

X glfk) ey:p{-ikx)E()])^^^^_^^y^.^{-hF^)dk 

°° (—nY 1°° 

+ (fc) exp(—ifccc) 

^(7l-72)(r-+l) + l^*( 7 .) 

X eyvp{-ikx)Eyy_^^^f^^^^.^^^{-bt'<^)dk 

+ / U*{k,t — f)exp{—ikx) 

r=0 -^0 J-OC 

>< (3.15) 

where 3 fJ( 7 i) > 0 , 5 ft( 72 ) > 0 and 5 R( 7 i — 72 ) > 0 , 6 = w + pip^^^k). 

4. A Particularly Interesting Case 

Theorem 2. Under the conditions of Theorem 1 with 1 < 71 < 2 replaced by 1 < < 1 and 

1 < 72 < 2 replaced 67 / 1 < 72 < 1 and following similar method, the solution of the following 
one-dimensional fraetional reaction-diffusion equation of fractional order 

D]^’^^N{x,t) + aDf^’^^N{x,t) =pxDe°‘N{x,t) -ujN{x,t) + U{x,t), (4.1) 

where uj,p,t > 0,x G R; a, 9, 71 , 72 , ^ 1 , <^2 are real parameters with the eonstraints 

1 < 7 i < 1,0 < Ji < 1; 1 < 72 < 2 = 1, 0 < (52 < 1, 0 < a < 2; (4.2) 

with the conditions 

^a- 5 i)(i- 7 i)^(^^ 0) = hi{x)-, N{x, 0) = h 2 {x)-, lim N{x, t) = 0. (4.3) 

|ai|—5-00 



is given by 



°° (—nV 

-/V(x, t) = ^ —-- / exp{—ikx) 

r =0 '^~°° 

, (-bt'^^)dk 

7l.7l+oi+(7l-72)i—7101' ' 

°° ( — nV t°° 

+ ^ —-- / exp(—ifcx) 

n 'J —oo 


^ -®7l.7l + (7l-72)r-+<52(l-72)( 

+ / C/*(fc,t-e)exp(-ifcx) 

r-=0 “^0 

where 9fJ(7i) > 0 , 5R(72) > 0 and 3fJ(7i — 72) > 0 , 6 = w + ri'ip^{k). 


(4.4) 


If we set 0 = 0, then by virtue of the relation (A14), the Riesz-Feller space derivative reduces 
to the Riesz derivative and Theorem 2 reduces to the following: 


Corollary 2.1. Under the eonditions of Theorem 2 with 9 = 0, the solution of the following 
fractional one-dimensional reaction-diffusion equation of fractional order 

t) + aD]^’^^N{x, t) = r],,Do°N{x, f) - ujN{x, t) + U{x, t), (4.4) 

with the conditions 

j(i-5i)(i-7i)^^^^ Q) _ /(1 “'52 )(i-72 ) _ /j 2 (x); lim N{x,t) = 0. (4.5) 

+ + |a:|—s-oo 

is given by 


°° (—nY 

Af(x,t) / t''^^-^^^'^hl{k)eyiY>{-ikx) 

t'o 

^ -®7lt7l+<5l+(7l-72)r-7l5l(“^^^')'^* 

°° (—nY 

+ ^ - / t^^^“'*'^^’'/i 2 (fc) exp(—ifecc) 

_n J — 00 


X ,(-6CMdfc 

7l.7l+(7l-72)i'+02(l-72) ' ' 

00 (_„\r ft poo 

+ y^- - - / ^^t+Yi-i2)r-i / u*{k,t-OeM-ikx) 

ZYf 27r Jo i -00 


r—Q 

-ir+l 


X _,,).(-6e)dMc, 

where lR( 7 i) > 0 , 5 ft( 72 ) > 0 and lR( 7 i — 72 ) > 0 , 6 = w + gif^^k). 


(4.6) 


Corollary 2.2. Under he conditions of Theorem 2 with uj = 0, the solution of the following 
one-dimensional fractional reaction-diffusion equation of fractional order 

Df^’^^N{x, t) + aD]^’^^N{x, t) = r],,Dg°‘N{x, t) + U{x, t), (4.7) 
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with the conditions 

4i”‘^')(i”^i)Af(x,0) =/ii(x); 4 i“'^")(i”^ 2 )^( 2 .^ 0 ) =/i 2 (a;); lim N{x,t)=t). (4.8) 

+ |a:|—s-oo 

is given by 


°° (—nY r°° 

N{x,t) = exp{—ikx) 

t'o 

°° (—nY f°° 

+ ^ —-- / exp(—ifccc) 

_n J — OO 


X A-b*f<^)dk 

7i.7i+(7i-72)»'+<52(1-72) ' ^ 


( — nY r 1'°° 

Y- - - / / U*{k,t-OeM-ikx) 

YY'n 27r Jo i-oo 


r=0 

Tir+l 


w/iere 3 fJ( 7 i) > 0 , 5 R( 72 ) > 0 and 5 R( 7 i — 72 ) > 0 , 6 * = r]^^{k). 


(4.9) 


If we set (5i = 0, ^2 = 0 then the generalized Riemann-Liouville fractional derivatives 
and reduce to Riemann-Liouville fractional derivatives and defined in 

(A5) and we arrive at the following: 

Corollary 2.3. Under the conditions of Theorem 2 with Si = 0 = 62 , the solution of one¬ 
dimensional reaction-diffusion equation of fractional order 

^^Dt^^N{x,t) + N{x,t) = ri^Dg°‘N{x,t) - u}N{x,t) + U{x,t), (4-10) 

for uj,r],t > 0,x € R; a, 0, 71 , 72 , are real parameters with the constraints with the conditions 
7j(7i-i)jy(2,^0)) = ei(x); 0+) = e2(x) 


lim N{x, t) = 0,t > 0. 

|ai|—)-oo 


(4.11) 


is given by 


°° (—nY 

lV(a;, t) = ^ —-- / e\(k) e-xp>{—ikx) 

r=o 


X 77;+,,. jdt 


nC'-'Y. 


r =0 


(z^ 

27r 


^(7i-72)r 


' —00 
-ir+l 


X el{k) eyip{-ikx)E:^+^^^^^^_^^y{-bU^)Ak 

( — nY U 

+ 'S^ - / ^"yi+Ui—yUr-i / u*{k,t — f)exp{—ikx) 

fr'n 27 r Jo J-00 


r—0 

pr+l 


X , (-br^)dkd£, 

7l.7l+(7l-72)»'' ^ 


(4.12) 
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where 3 fJ( 7 i) > 0 , 3 ft(j 2 ) > 0 and 5 R( 7 i — 72 ) > 0 , 6 = w + 

If (5i = ^2 = 1, the generalized Riemann-Liouville fractional derivatives and 

reduce respectively to the Caputo fractional derivatives 0 Dt'^^ and 0 defined by (A6) and 
we arrive at the following: 


Corollary 2.4. Under he conditions of Theorem 2 with hi{x) = h 2 {x) = h(x) and di = 62 = 1, 
the solution of the following one-dimensional fractional reaction-diffusion equation of fractional 
order 


^Et'''N(x,t) + 

Oq Dt'^^N{x,t) = 'nxT)e°‘N{x,t) — ujN{x,t) + U{x,t) 

(4.13) 

with the constraints 

N{x, 0+) = h{x), lim N{x, t) = 0. 

(4.14) 


|ai|—>-cxD 


is given by the following: 




°° (—aY C°° 

N{x, t) = I t^^^~'*^^'^h*{k)exp{—ikx) 


r=0 

7r+l 


xE^+j , ^J-bE^)dk 

7i.(7i-72))'+1' ' 




(-a)^ 


r =0 J-00 


^(7i-72)(r-+l) 


X h*{k) 

“ ( — aY U 

+ / U*{k,t — f)exp{—ikx) 

7^0 "'o 

xE''+^^, , (-bC^)dkd£, 

7l.7l + (7l-72)»'' ^ ’’ 


(4.15) 


where 9fJ(7i) > 0,5R(72) > 0 and lR( 7 i — 72 ) > 0, 6 = w + ri'ipf^{k). 

For w = 0 (4.13) reduces to one given by Saxena et al. [50]. If we further set h{x) = (5(x), 
where S(x) is the Dirac-delta function, then Theorem 2 can be written in the following form: 


Corollary 2.5. Under the conditions of Theorem 2, the fundamental solution of the following 
one-dimensional reaction-diffusion equation of fractional order 

Q N{x,t) + Oq N{x,t) = 'ilxDe°‘N{x,t) —u)N{x,t) + U{x,t) (4-16) 

with the constraints 

N(x, 0+) = d(x), lim A^(a;, t) = 0, t > 0. (4-17) 

|ai|—)-cxD 
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is given by the following: 


°° (—nY 

N{x,t) = '^^— - / exp(—ifca;) 


r =0 


{-ay 

27r 


^(7i-72)(r-+l) 


' —oo 
^r +1 


“ (—nY Y r°° 

+ / U*{k,t — ^)exp{—ikx) 

7^0 "'o 

x£;’'+^^, , (-6P)dfcd7 

7l.7l+(7l-72)»'' ’ 


(4.18) 


where 9 fJ( 7 i) > 0 , 5 R( 72 ) > 0 and 5 R( 7 i — 72 ) > 0,6 = w + rjip^^k). 

For w = 0, (4.17) yields another result given by Saxena et al. [50]. As a concluding remark, 
it is interesting to observe that Theorem 2 also hold true if instead of one Riesz-Feller derivative, 
we consider a finite number of Feller derivatives. This result is given in Theorem 3 in the next 
section. 


5. Several Riesz-Feller Space Fractional Derivatives 

Theorem 3. Consider the one-dimensional fractional reaction-diffusion equation of fractional 
order 

m 

oy’^^Nix, t) + aD^’^^Nix, t) = N{x, t)] - cjN{x, t) + U{x, t), (5.1) 

1=1 

where uj,r]j,t > € R; aj,0j,j = 1,..., m, 71 , 72 , Si, 62 are real parameters with the constraints 

1 < 7 i < 2,0 < (5i < 1; 1 < 72 < 2,0 < (^2 < 1,0 < Oj < 2, j = 1, ...,to; (5.2) 

andD^’^^ are the generalized Riemann-Liouville fractional derivative operators defined 
by (A9) with the conditions 

4E')("-7^)^7,0) = /2(X); ^ 

lim N{x,t) = 0. (5.3) 

|ai|—)-oo 

Further, u is a constant with reaction term, xDgC^ are the Riesz-Feller fractional derivatives of 
order aj and symmetry 9j defined by (All) with |6*j| < min(Q;j,2 — aj),) = 1, ...,m,r]j > 0, j = 
l,...,m are the diffusion constants and U{x,f) is a nonlinear function belonging to the area of 
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reaction-diffusion. Then the solution of (5.1), under the above conditions, is given by 


°° (—oY 


^ -®7lt7l + (7l-72)r-+<5i(2-7i)-l( 

_|_ ^ 7 i +<^ i ( 2 - 7 i )-1 ' (~'^) ^(71-72)1’ 


r—0 


27r 


°° (—nY 

+ ^ — - / (A:) exp(—ifcx) 

_^q^7i+-52(2 - 72)-1 ^ (Z^ f t(7i-72)r-^*(^) 

r=0 ^ d-oo 

X exp(—^ A—bt^^)dk 

( — flY Y f°° 

+ - / ^ 7 i+( 7 i- 72 )»'-i / U*{k,t — ()exp{—ikx) 

t^o ^TT Jo J-00 

x£;’'+^^, , (-6^'^MdfcdC, 

7l.7l + (7l-72)r' 7 ; S) 


w/iere K( 7 i) > 0 , 5 R( 72 ) > 0 and ^(71 - 72 ) > 0, 6 = w + J^jLi Vj^’a) {k). 

6. Special Cases of Theorem 3 


(5.4) 


Corollary 3.1. Under the conditions of Theorem 3, with u) = 0, the one-dimensional fractional 
reaction-diffusion equation of fractional order 

m 

D(^’^^N{x,t)+aDf^’^^N{x,t) = N{x,t)] + U{x,t), (6.1) 

i=i 

where rij,t > 0,x € R;aj,0j,j = 1,..., m, 71 , 72 , di, (52 are real parameters with the constraints 

1 < 7 i < 2 , 0 < di < 1 ; 1 < 72 < 2,0 < (^2 < 1,0 < Oj < 2 , J = 1 ,..., m; ( 6 . 2 ) 

are the generalized Riemann-Liouville fractional derivative operators defined 
by (A9) with the conditions 

4r'')(2-7^)iV(x,0) = hix); Aja-^i)(2-7d^(^^0+) = 5 i(x) 

4r'^)('-72)lV(x,0) = hix); = g^ix), 

lim N{x,t) = 0. (6.3) 

|ai|—)-oo 

Further, xDgT^ are the Riesz-Feller fractional derivatives of order aj and symmetry 9j \6j\ < 
min(Qfj,2 — aj),j = l,...,m,gj > 0,j = l,...,m are the diffusion constants and U{x,t) is a 
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nonlinear function belonging to the area of reaction-diffusion. Then the solution of (6.1), under 
the above conditions, is given by 


OO /_ \r poo 

N{x,t) = exp(—ifcx) 


_j_ ^ 7 i + i 5 i ( 2 —71) — 1 ^^ (__^)_ ^(71-72)1’ 


r—0 


27T 


°° (—nY 1 '°° 

+ ^ —-- / exp(—ifcx) 

r=o 

xE^+^^, , , J-W^)dk 

7i.7i+(7i-72)»'+02(2-72)-1' ' 

00 / \J. «00 

_^Q^7l+-52(2 - 72)-iy^ / t^'^^-'f^>g*{k) 

f-i 27r 


°° (—nY /■* 1°° 

+ Y- - - / / U*{k,t-f)eM-ikx) 

ZY( 27r Jo i -00 


r=0 

T-r+l 




(6.4) 


w/iere 9fJ(7i) > 0,5R(72) > 0 and K( 7 i - 72 ) > 0, 6 = Vj''P% (k)- 

When 9i = 0 = 02 = ... = 9m the Riesz-Feller space derivatives reduce to Riesz Fractional 
derivative, then by virtue of the relation (A14), there holds the following: 


Corollary 3.2. Under the conditions of Theorem 3, the fractional reaction-diffusion equation 
of fractional order 


Df^’^^N{x, t) + aDf^’^^N{x, t) = YyijxDo°‘^N{x, t)] - ojN{x, t) + U{x, t), (6.5) 

j=i 

where uj,t > 0,x € Tl;r]j > 0,aj,j = 1,..., m, 71 , 72 , (52 are real parameters with the con¬ 

straints 

1 < 7 i < 2,0 < < 1; 1 < 72 < 2,0 < (^2 < 1,0 < Oj < 2, j = 1, ...,m; (6.6) 

andD]-^’^^ are the generalized Riemann-Liouville fractional derivative operators defined 
by (A9) with the conditions 

4+”'^')(2”7')lV(a;,0) = fi{x); N{x,0+) = gi^x) 

4r'^)^2-72)^(^^0) = / 2 (x); = g 2 {x), 

lim A^(a;,t) = 0. (6.7) 

|ai| —5-00 
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Then the solution of (6.1), under the above conditions, is given by 


oo (_„Y roo 

N{x,t) = exp(—ifcx) 

r=0 ^ d-oo 

^ -®7lt7l+(7l-72)r-+<5i(2-7i)-l(“^*^^')d* 

_|_ ^ 7 i +<^ i ( 2 - 7 i )-1 (~'^) f ^{'n-'t2)r 

_n ^TT 


°° (—nY f°° 

_^Q^7i+52(2-72)-2y-i^ / t(^i-'^=)’'/2(fc)exp(-ifcx) 

t'o 

^ ^ J-b*t'^^)dk 

7i57i + (7i—72)r*+<32(2—72) —1 ^ ' 

(—nY 

_^Q^7i+-52(2 - 72)-iy^ -L / g*{k) 

Yo 

°° ('—o'l’' Y r°° 

y / p + (7i-72)r-l / _^)exp(-ifcx) 

_n -/o J-OO 


r =0 

T-r+l 




( 6 . 8 ) 


w/iere K( 7 i) > 0,5R(72) > 0 and 5 R( 7 i - 72 ) > 0 , 6 * = w + X)^i VjY) (^)- 

If we set Si = 62 = 0, then the generalized Riemann-Liouville fractional derivatives 
and defined by (A5) and we arrive at the following result: 


Corollary 3.3. Consider the one-dimensional fractional reaction-diffusion equation of frac¬ 
tional order 

m 

YdY'"N{x, t) + aYDY'"N{x, t) = YYo-Dep N{x, t)] - ojN{x, t) + U{x, t), (6.9) 

i=i 

where uj,t > 0,x € R;r]j > 0 ,aj,j = 1,...,m, 71,72 are real parameters with the constraints 

1 < 7i < 2,1 < 72 < 2,0 < Oj < 2, j = 1,..., m; (6.10) 

and ^^6 generalized Riemann-Liouville fractional derivative operators 

defined by (A5) with the conditions 

Df^-^N{x,0) = fi{xy,D(^-^N{x,0+) = gi{x) 

D^-^Nix, 0 ) = / 2 (x); D^-^Nix, 0 ) = g 2 {x), 
lim N{x,t)=0 (6.11) 

|ai|—^00 


where D^jj — 2)N{x,0) denotes the (jj — 2)th derivative of N{x,t) evaluated at t = 0, to is a 
constant with reaction terms, xDgC^ are the Riesz-Feller space fractional derivatives f order aj 
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and symmetries Oj; \0j\ < min(Q;j,2 — aj),j = defined by (All), ijj > 0 are dijfusion 

constants and U(x, t) is a conlinear function belonging to the area of reaction-diffusion. Then 
the solution of (6.9), under the above conditions, is given by 


N{x,t)=t^^-^^ 


r—0 


i-ay 

27T 


^(71 exp(—ifcx) 




., 7 .- 1 ^ 


(-a)" /■“ 


r=0 


27r 


^(7l-72)r 


' —00 
Tir+l 


X g({k) exp{-ikx)E(^+^^^^^^_^^^^{-qt^^)dk 




r—0 


i-aY 

27r / _ 


^(71 


X , J-qEY<^k 

7l.7l + (7l-72)»'-l' ^ ' 




r—0 


i-ay 

27T 


i(7l-72)r^*(^) 


X exp(—ikx)E^'^ . {—qE^)dk 

t'y ! 71)71 + (71—72)r*'' ^ ’ 

( — nV 

+ —-— / / U'^(^k,t — ^)exp{—ikx) 

_n ^TT Jq J-oo 


r—0 

pr+1 




( 6 . 12 ) 


where 3 fJ( 7 i) > 0 , 5 ?( 72 ) > 0 and 5 R( 7 i — 72 ) > 0 ,q = ^ 2^=1 • 

For w = 0 the above result reduces to the result given by Saxena et al. [48]. If we set 
di = d 2 = 1 then the generalized Riemann-Liouville fractional derivatives and 

reduce respectively to the Caputo fractional derivatives 0 dY^ and (f dY^ defined by (A 6 ) and 
we arrive at the following: 


Corollary 3.4. Under the conditions of Theorem 3, the solution of one-dimensional fractional 
reaction-diffusion equation of fractional order 


((DY^N{x,t) + Oq dY^N{ x,t) = ['^rijxD 0 YyA[{x,t) - ujN{x,t) + 17(x,t) (6.13) 

j=i 

where uj,t > 0; rjj > 0, aj, 0 j,j = 1,..., m, 71 , 72 are real parameters with the conditions 1 < 71 < 
2,1 < 72 < 2,0 < Qfj <2 and with the constraints 

N{x, 0+) = f{x); -^Nix, 0+) = g{x) (6.14) 

dx 

and other quantities are as defined in Theorem 3. Then the solution of (6.13) is given by the 
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following: 


°° (—oY 

N{x,t) = '^—^ - / f*{k) exp{—ikx) 


r =0 

^r +1 


i-ar 

27T 


fhi exp(—ifcx) 


°° (—nY f°° 

+ a^—-- / f*{k)exp{—ikx) 


^E 

r—O 


i-ay 

27T 


^(7i-72)(r-+l) + l 


“ ( — aY Y 

+ / U*{k,t — f)e^pi—ikx) 

ZY'n 27r Jo J-oo 


r =0 

7ir+l 


X _,,)7-67'^7dM7 


w/iere 3 fJ( 7 i) > 0 , 5 i( 72 ) > 0 and 5 R( 7 i - 72 ) > 0 , 6 * = w + (k). 

7 . Conclusions 


(6.15) 


In this paper, the authors have presented the solutions of three one-dimensional fractional 
reaction-diffusion equations of fractional orders associated with generalized Riemann-Liouville 
derivative as the time derivative due to Hilfer et al. [23] and Riesz-Feller derivatives as the space 
derivatives. The results obtained provide an elegant extension of the solution of one-dimensional 
fractional reaction-diffusion equations associated with Caputo fractional derivatives as the time 
derivatives and Riesz-Feller fractional derivatives as the space derivatives [50]. The results are 
obtained in terms of generalized Mittag-Leffler functions. Further, the derived results include, 
as a special case, the results for the solution of space-time fractional reaction-diffusion systems 
associated with a generalized Riemann-Liouville fractional derivative given by the authors [49], 
which itself is a generalization of the fundamental solution of space-time fractional diffusion 
given by Mainardi et al. [29]. The solution of the equations (2.1),(4.1), and (5.1) are obtained 
in terms of generalized Mittag-Leffler functions. The importance of the results obtained in this 
paper further lies in the fact that due to presence of the modified Hilfer fractional derivative 
[23], results for Riemann-Liouville fractional derivative and Caputo fractional derivatives can 
be deduced as special cases by taking n = 0 and (5i = 0 and i 52 = 1 respectively. 
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Appendix: Mathematical Preliminaries 

A generalization of the Mittag-Lefher function (Mittag-Lefher [32,33]) 

oo „ 

Ec.iz) = J2—^^,aGC,^{a)>0 (Al) 

1 (na -I- 1) 

n=0 

was introduced by Wiman [64] in the generalized form 

OO 

> 0- (^2) 

A further generalization of the Mittag-Lefher function is given by Prabhakar [39] in the form 

Tl 

ElA^) = E r/ > 0, (A3) 

n—0 ^ ' 

where the Pochhammer symbol is given by 

{a)n = a{a -I- 1)...(a -I- n — 1), a 0, (a)o = 1. 

The main results of the Mittag-Lefher functions dehned by (Al) and (A2) are available in the 
handbook of Erdelyi et al. [9, Section 18.1] and the monographs written by Dzherbashyan [5,6]. 
The left-sided Riemann-Liouville fractional integral of order v is dehned by Samko et al. [40], 
Kilbas et al. [26], Mathai et al. [31] as 

ADt ''N{x,t) = -—— / {t — u)''~^N{x,u)du,t > OAA > ^- (^4) 

r(f^) Jo 


,na. 
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The left-sided Riemann-Liouville fractional derivative of order a is defined as 

“A^iV(x, t) = ( >0,n = [»(/x)] + 1, (A5) 

where [x] represents the greatest integer in the real number x. Caputo fractional derivative 
operator (Capuato[l]) is defined in the form 


C 

0 




1 /•* /(’")(x,T)dT 

r(TO — a) Jq (t — 

—--, mboxija = m 


, m — 1 < a < m3?(a) > 0, m € 


(AG) 

{A7) 


where ^^f(x,t) is the mth partial derivative of f(x,t) with respect to t. When there is no 
confusion, then the Caputo operator CoDt°‘ will be simply denoted by oDf. A generalization of 
the Riemann-Liouville fractional derivative operator (A5) as well as Caputo fractional derivative 
operator (A6) is given by Hilfer [22] by introducing a left-sided fractional derivative operator of 
two parameters of order 0 < /i < 1 and type 0 < 1 in the form 


Dil':N{x,t) 



m 


For 1 ^ = 0, (A8) reduces to the classical Riemann-Liouville fractional derivative operator (A4). 
On the other hand, for = 1 it yields the Caputo fractional derivative operator defined by (A5). 


Note Al: The derivative defined by (A8) also occurs in recent papers by Hilfer et al. [23], 
Srivastava et al. [59], Tomovski [60], Tomovski [61,62] and Saxena et al. [52]. Recently, Hilfer 
operator dehned by (A8) is rewritten in a more general form (Hilfer et al. [23]) as 


D^a:N{x,t) 


An 


(A9) 


Where n — l</r<n, 0<^<1. The Laplace transform of the above operator (A9) is given by 
Tomovski [59,p.l7] in the following form: 


n—1 

L[D^+N{x, ty,s]= sf^N{x, s)-J2 

k^O 

(AlO) 

where n — 1 < ^ < n,n € N,0 < i/ < 1. Following Feller [10,11], it is conventional to define 
the Riesz-Feller space fractional derivative of order a and skewness 9 in terms of its Fourier 
transform as 




(All) 
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where 


^a(^) = l^r exp[i(si 5 nA:)—], 0 < a <2,\9\ < min(a, 2 — a). (^12) 

When 0 = 0, we have a symmetric operator with respect to x, that can be interpreted as 

-V = (-<13) 

This can be formally deduced by writing —(fc)“ = —For 0 = 0, we also have 

F{,,Do‘^f(x);k} = -l^rr (fc). (AU) 

For 0 < a < 2 and |0| < min(a, 2 — a), the Riesz-Feller derivative can be shown to possess the 
following integral representation in x domain 


,De‘^f{x) = Iii±^[sin((a + 0)^/2) H 


+ sin((Q; — 0 ) 7 t / 2 ) ( 
Jo 


JO 

f{x -0- f{x) 


^1+0 


-dC]. 


For 0 = 0, the Riesz-Feller fractional derivative becomes the Riesz fractional derivative of order 
a for 1 < a < 2 defined by analytic continuation in the whole range 0<Q;<2,ay^l (see 
Gorenflo and Mainardi [14]) as 

= -A[/+“ - /!“] (A15) 


where 


A = 


1 


2cos(Q;7r/2) 


I -'it 


d2 

fx2 


(^16) 


The Weyl fractional integral operators are defined in the monograph by Samko et al. [40] as 


iI^N)ix) = ^ I Jx - > 0 

1 

{I^-N){x) = y (^ - P>o. (4117) 


Note A2; We note that xDq^ is a pseudo differential operator. In particular, we have 


= ITT’ ^ 

ax^ ax 

The following result given by Saxena et al. [44] is also required: 


(A18) 


l-M- 


op-1 


as^ 


-;t] = t° 


r=0 


A / a,Q:+(a —p)r—p+1 A / 


(A19) 


where 5ft(a) > 0, K(/3) > 0, 1 R(q; — /3) > 0,5ft(a — p) > 0,5R(s) > 0, < 1 and E'^ ^(z) is the 

generalized Mittag-Leffler function defined by (A3). 
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